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Abstract 

By explicitly allowing for topology to change as a function of time, two-dimensional 
quantum gravity defined through causal dynamical triangulations gives rise to a 
new continuum string field theory. Within a matrix-model formulation we show 
that - rather remarkably - the associated sum over all genera can be performed in 
closed form, leading to a nonperturbative definition of CDT string field theory. We 
also obtain explicit formulas for the n-loop correlation functions. Our construction 
exhibits interesting parallels with previous, purely Euclidean treatments. 
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1 Introduction 



By general acknowledgement, string theory is in need of a genuinely nonperturba- 
tive definition if it is to serve as a fundamental theory of all interactions. String 
theory's perturbative expansion does not serve this purpose, since it is only asymp- 
totic and becomes useless when the string coupling approaches order one. It is thus 
of great interest to study toy models like non-critical string theory for which one 
can define the sum over all genera, a possibility first outlined in the seminal papers 
[Tl[2],[3]. Unfortunately, the results of such summations are plagued by ambiguities: 
even in the simplest cases, like pure two-dimensional Euclidean quantum gravity 
(c = non-critical string theory) and all unitary non-critical string theories, the 
"nonperturbative" partition function is not Borel-summable [H El EJ [7] (see also [8] 
for a recent lucid discussion). Attempts to perform the summation in these theories 
typically result in complex partition functions or, perhaps more accurately, con- 
structions which do obtain a real partition function are often guided by the desire 
for a real outcome, rather than by applying stringent physical criteria. (For other 
approaches see, for example, [9|). 

Given this somewhat unsatisfactory state of affairs, any new physical model for 
which one has analytic control over the sum over topologies is potentially valuable, 
to understand possible ambiguities in the summation and their genericity, to explore 
their physical meaning, and to see whether new insights can be gained into the nature 
of surfaces of infinite genus, a topic which has maybe not received the attention it 
deserves. 

The present article deals with just such a new system in two spacetime dimen- 
sions. Its origin lies in the attempt to formulate a nonperturbative theory of four- 
dimensional quantum gravity based on conventional quantum field theory, more 
precisely, the path integral approach, in the framework of "Causal Dynamical Tri- 
angulations" (or "CDT" for short), which has led to a variety of intriguing and 
unprecedented results [HI [T2]. However, the construction can in principle be ap- 
plied in any dimension, leading to nonperturbative models of dynamical quantum 
geometry with causal, Lorentzian properties. The case of two spacetime dimensions 
is a particularly attractive testing ground because of the availability of a whole range 
of analytical tools. The original, strictly causal and purely gravitational model was 
solved exactly in 2d [13], and has since been generalized by the inclusion of matter 
[m [15] and isolated (cap or branching) points where causahty is violated. The latter 
has culminated in the recent definition of a fully fledged CDT string field theory (in 
zero-dimensional target space), which is the subject of the remainder of the present 
work. 

In the mid-nineties, Kawai, Ishibashi and collaborators developed a string field 
theory for non-critical strings [161 13 UHl [HI |20l |2T]. While its original starting 
point was the explicit realization of non-critical string theory in terms of dynamical 
triangulations (of the Euclidean variety), it eventually was formulated entirely in the 
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language of continuum field theory. In [22], [23] we repeated this continuum analysis 
for the two-dimensional CDT model of quantum gravity. We demonstrated how 
it is connected to the conventional matrix model of 2d Euclidean quantum gravity 
[211 125] and how the continuum Dyson-Schwinger equations of the CDT string field 
theory are related to a particular matrix modelQ In what follows, we will use this 
matrix-model formulation to perform the sum over all topologies in the continuum 
CDT string field theory and compute a number of associated geometric observables, 
so-called loop-loop correlation functions. It turns out that this can be done in a 
surprisingly simple manner. 

The remainder of this article is organized as follows: in section [2] we introduce the 
CDT matrix model and related observables, in section [3] we show how to calculate 
these quantities summed over all genera, in section H] we put our results into the 
context of previous Euclidean models, and finally, in section |5] we discuss the nature 
of our nonperturbative results. 

2 The CDT matrix model 

The approach of causal dynamical triangulations aims to define quantum gravity 
nonperturbatively via a path integral over geometries [Qhv] , 

Z(Gn,A) = j e-^[^-l, (1) 

where 5* denotes the (Euclidean) Einstein-Hilbert action of gravity, Gn is Newton's 
constant and A the cosmological constant. We are presently interested in the case 
of two spacetime dimensions for which the action is given by 

S[g^,u\ = -^^^ / y^detg^, R + \j(fi ^d^t^. (2) 

One proceeds in two steps, first setting up the CDT lattice regularization of the path 
integral in Lorentzian signature, and then rotating each individual triangulation 
in the Lorentzian ensemble to Euclidean signature. The resulting Euclidean path 
integral (jl]) differs from a standard Euclidean path integral since by construction 
each geometric "path" has a time-foliation and thus a "memory" of the causal 
properties it possessed before the Wick rotation. We refer the reader to the original 
article [IB] for technical details. 

The observables we will discuss in this article are generalized Hartle-Hawking 
amplitudes involving 2d geometries with a fixed number of boundaries. The sim- 
plest one is that for a single boundary, with all boundary points sharing a common 

"'^Since no scaling limit is needed in this matrix model, the situation is similar to what happens 
in topological 2d gravity where the Kontsevich matrix model directly describes the continuum 
gravity theory [26l[27l[28] (see [29] for a recent description in terms of D-branes). As will become 
clear below, this similarity is not accidental. 
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time (in the time- foliation referred to above). This so-called Hartle-Hawking or disc 
amplitude W{1) is the amplitude that the spatial boundary (consisting of a closed, 
one-dimensional loop) has length /. Depending on whether one thinks of the bound- 
ary as being initial or final, it is the amplitude of the universe of length / vanishing 
into or coming from "nothing" . A generalized such amplitude is one where the spa- 
tial boundary (still at equal time) consists of n disconnected loop components with 
specified lengths li, . . . ,ln- Spacetime may remain disconnected or reconnect in var- 
ious ways at other times. We denote the corresponding amplitudes Wd{h, ■ ■ ■ Jn), 
where the subscript "d" indicates that for n > 1 it includes the possibility of having 
disconnected universes. The connected component of Wdih, ■ ■ ■ ,ln) "will be called 
Wc{h, ■ ■ ■ ,ln)- To lowest order n = 1 we have Wd{l) = Wc{l) = W{1), while for n = 2 

W,{h, h) = WSi. h) - W{h)W{h), (3) 

which generalizes to higher-loop correlators in standard fashion. 

Instead of considering a situation where the boundary lengths li are kept fixed, 
one can introduce by hand boundary cosmological constants Xj and include corre- 
sponding terms Xi ■ U in the action. Now the lengths U are allowed to fluctuate, but 
their average values are determined by the values of the Xi. This leads to another 
kind of generalized Hartle-Hawking wave functions, denoted by Wd{xi, . . . ,a;„) and 
Wc{xi, . . . , Xn), and related to the previous VFc,d(^i, • • • , ^n) by Laplace transforma- 
tion according to 



/•oo 

W,,d{xi, . . . , x„) = / d/i ■ ■ ■ d/„ ^^^^^^^^ _ _ _ ^ 

Jo 



(4) 



The Dyson-Schwinger equations, an infinite set of coupled equations relating 
Wc{xi, . . . ,Xn) for different n, follow from the continuum CDT string field theory 
as developed in [22]. They are obtained from the cubic matrix model model 



ZNi9,\) = jdM exp 



-— tr (\M~- 
9 \ 3 



(5) 



where M is an x Hermitian matrix. In this matrix model the observables 
introduced above take the form 









^tr 












^tr 







= ^(t.l^V..tr(-l^l) , (7) 



connected 



where the expectation value of an operator O is defined by 



(0(M)) = ^ } f dM OiM) exp -— tr (AM - - MM 

ZN{g,X)J L ^ V 3 y 
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Of course an expression like ([5]) is formal and should always be understood as 
a suitable expansion in powers of M . This is most clearly exhibited by making the 
variable change 

M =Vt\Y ^=^ = e-l/^^ (9) 



where t is a dimensionless coupling constant, which is related to the coupling g of 
eq. ([5]) and the gravitational coupling Gn as specified. In terms of the matrix Y and 
the coupling t the matrix integral reads 



Z^(t(^?,A),A) = exp(^') (tA)^^/2 



dY exp 



-iVtr Y 



Y' 



(10) 



which after expanding exp{\^Y^ /3) and performing the Gaussian integrals becomes 
a formal power expansion in t. The coefficient of in this power series is positive, in 
agreement with the interpretation that it represents the number of 2d surfaces with 
a time foliation with n degenerate points (see p2| for details). Let us emphasize 
that the pre-f actor, which does not have a power expansion in t, will cancel in any 
expectation values of observables and should not really be considered part of the 
partition function. In fact we could have defined the partition function as 



ZNit) 



J dY exp -A^tr (y^ - ^Y^^ 



f dY exp [-iVtr Y^] 



(11) 



which would then only contain positive powers of t. However, we find it convenient 
to use ffTOj) in the rest of this paper, with the understanding that the exponential 
growth of Z]^{t) for t ^ should not be considered as a refiection of unphysical 
behaviour, but is merely a consequence of a particular choice of normalization. 



3 Summing over all genera 

In reference [21] the matrix model ^ was related to the CDT Dyson- Schwinger 
equations by (i) introducing into the latter an expansion parameter a, which kept 
track of the genus of the two-dimensional spacetime, and (ii) identifying this param- 
eter with l/N'^i where N is the size of the matrix in the matrix integral. The 
expansion of our matrix model therefore plays a role similar to the 1/A^-expansion 
originally introduced by 't Hooft [30j: it reorganizes an asymptotic expansion in 
a coupling constant (t in our case) into convergent sub-summations in which the 
fcth summand appears with a coefficient N~'^'^ . In QCD applications, the physically 
relevant value is = 3, to which the leading-order terms in the large iV-expansion 
can under favourable circumstances give a reasonable approximation. 

As we will see, for the purposes of solving our string field-theoretic model non- 
perturbatively, an additional expansion in inverse powers of (and thus an identi- 
fication of the contributions at each particular genus) is neither essential nor does 
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it provide any new insights. This means that we will consider the entire sum over 
topologies "in one go", which simply amounts to setting = 1, upon which the 
matrix integral (jSD reduces to the ordinary integral 



dm exp 



Am 



m 



(14) 



while the observables (ED can be written as 



Z{gA) 



dm 



exp 



- (Am — I m^) 



{xi — m) ■ ■ ■ {xn — m) 



(15) 



Again, these integrals should be understood as formal power series in the dimen- 
sionless variable t as mentioned below eq. ffTOl) . Any choice of an integration contour 
which makes the integral well defined and reproduces the formal power series is a 
potential nonperturbative definition of these observables. However, different con- 
tours might produce different nonperturbative contributions (i.e. which cannot be 
expanded in powers of t), and there may even be nonperturbative contributions 
which are not captured by any choice of integration contour. As usual in such 
situations, additional physics input is needed to fix these contributions. 

To illustrate the point, let us start by evaluating the partition function given in 
( IT^ . We have to decide on an integration path in the complex plane in order to 
define the integral. One possibility is to take a path along the negative axis and 



^Starting from a matrix integral for N x A'^- matrices like (O, performing a formal expansion in 
(matrix) powers commutes with setting N — 1, as follows from the following property of expectation 
values of products of traces, which holds for any n — 1,2,3,... and any set of non-negative integers 
{nfc}, = 1, . . . , 2n, such that X]fc=i ^fc — 2ri. For any particular choice of such numbers, consider 



2« ^ -I , f dM e- HfelLi f trM"^ /N 
(TT ( _LtrM"M) = ^-^ 



E 



(12) 



where the last equation defines the numbers ujr, 
expectation value. Now, we have that 



as coefficients in the power expansion in N of the 



E 



Lo„, ^ (2n-l)!! 



(13) 



independent of the choice of partition {rifc}. The number (2n — 1)!! simply counts the "Wick 



contractions" of x which we could have obtained directly as the expectation value 



evaluated 



with a one-dimensional Gaussian measure. In the model at hand, we will calculate sums of the 
form X]m=-n '^"i directly, since we are summing over all genera without introducing an additional 
coupling constant for the genus expansion. In other words, the dimensionless coupling constant t in 
this case already contains the information about the splitting and joining of the surfaces, and the 
coefficient of t^ contains contributions from 2d geometries whose genus ranges between and [fc/2]. 
We cannot disentangle these contributions further unless we introduce N as an extra parameter. 
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then along either the positive or the negative imaginary axis. The corresponding 
integrals are 

Z{g, A) = t'/^F±{t-^/^), F±(r2/3) = 2n e^'^/^^Aiir^/'e^^^'/^), (16) 

where Ai denotes the Airy function. Both F± have the same asymptotic expansion 
in t, with positive coefficients. Had we chosen the integration path entirely along 
the imaginary axis we would have obtained (27ri times) Ai(t~^/^), but this has an 
asymptotic expansion in t with coefficients of oscillating sign, which is at odds with 
its interpretation as a probability amplitude. In the notation of ^31j we have 

Fi{z) = 7r(Bi{z)±tAi{z)^, (17) 

from which one deduces immediately that the functions F-j-(t~^/^) are not real. How- 
ever, since Bi(t~^/^) grows like est for small t while Ai(t^^/'^) falls off like e~3t , their 
imaginary parts are exponentially small in 1/t compared to the real part, and there- 
fore do not contribute to the asymptotic expansion in t. An obvious way to define 
a partition function which is real and shares the same asymptotic expansion is by 
symmetrization, 

l(F+ + F_)=7rBi. (18) 

The situation parallels the one encountered in the double scaling limit of the "old" 
matrix model [6], and discussed in detail in [8], but is less complicated. We will 
return to a discussion of this in the next section. 

Presently, let us collectively denote by F{z) any of the functions F±{z) or 7rBi(2;), 
leading to the tentative identification 

Z{g,X) = VXt^^^F(t-^/'''^, F"{z)=zF{z), (19) 

where we have included the differential equation satisfied by the Airy functions for 
later reference. In preparation for the computation of the observables Waixi, . . . , a;„) 
we introduce the dimensionless variables 

x = y^/X, m = g^/^(3, Wd{xi, . . . ,Xn) = X-'^^^Mvi, ■ ■ ■ ^Vn)- (20) 
Assuming j/^ > 0, we can write 

1 

-—^ = da exp [- {y - t'/'p) a] . (21) 

We can use this identity to re-express the pole terms in eq. f|T5l) to obtain the integral 
representation 

My,, ...,yn)= l[da, e-(^^-^+-+^"°") ^ 

Jo t\ F t-^/^^ 
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for the amplitude with dimensionless arguments. From the exphcit expression of 
the Laplace transform, eq. (jl]), we can now read off the generalized Hartle-Hawking 
amplitude as function of the boundary lengths, 

(23) 

For the special case n = 1 we find 

for the disc amplitude, together with the remarkable relatiorH 

W,{h,...,Q = W{h + --- + Q. (25) 
By Laplace transformation this formula implies the relation 

W,{xu . . . , x„) = n-T""! V (26) 

Before turning to a discussion of the nonperturbative expression for W{1) we have 
just derived, let us remark that the asymptotic expansion in t of course agrees with 
that obtained by recursively solving the CDT Dyson- Schwinger equations. Using 
the standard asymptotic expansion of the Airy function [31] one obtains 



W{1) = e-v^' e*'^^*'^') ^^=''''1S1~^'^^^'"'~\ (27) 

where the coefficients Ck are given by cq = 1, Ck = ^ (f)'^ (|)^ (1)^' ^ > 0- (123), 
we have rearranged the exponential factors to exhibit the exponential fall-off in the 
length variable I, multiplied by a term containing the function 



{2i 



which has an expansion in positive powers of t. 

Finally, let us derive an expression for the amplitude W{x). Since F[z) satisfies 
the Airy differential equation, we have the identity 



a). (29) 



■^Formula (|25p has a structure quite similar to the one encountered in the "old" matrix model 
[321133] . where again the multi-loop correlator is "almost" a function of + • • • + /„ only. 
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Inserting this into eq. (122|) and performing the partial integrations, we obtain a 
first-order differential equation for w{y), namely, 



F U-2/3 



(30) 



Its solution is given by 



dv e 



{v-l-f{t)), 



(31) 



where the function f(t) is defined by 



m = — - 1. 



(32) 



Performing a variable shift f = ?/ + in the integral in (!3T|) . w(?/) is conveniently 
written as 



wiy) = d^e 



{y-l)-f(t)+t^ 



(33) 



From (133|) it follows by expanding the exponential containing t that w{y) has an 
asymptotic expansion in t (for ?/>!). The same expansion represents w{y) by 
an expansion in inverse powers of {y + 1), corresponding to the expansion fl271) . 
Explicitly, one finds 



w{y) 



1 y + 3 ^/ 9x 



y + 1 % + l)3 



(34) 



which, as already stated, of course agrees with the perturbative expansion of w{y) 
derived previously in CDT string field theory [22j . 



4 Relation with other models 

We should perhaps not be too surprised to meet the Airy function as part of our 
nonperturbative analysis, since it has already appeared previously in non-critical 
string theory, more specifically, in the so-called 2d topological quantum gravity. 
This theory can be described in two ways. On the one hand, a set of observables of 
the theory, the intersection indices of Riemann surfaces, can be calculated by the 
Kontsevich matrix integral, which is the matrix generalization of the Airy function. 
This is a cubic integral like ([5]), but with A a matrix, which gives it more structure 
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than ([5]) and allows for the calculation of the intersection indices. Analogous to our 
case dSD, it is a matrix representation of a continuum theory, namely, 2d topological 
quantum gravity. 

On the other hand, 2d topological quantum gravity also has a "conventional" 
one-matrix representation, which in fact is the simplest one possible. Recall that 
the [p, q) minimal conformal field theories coupled to 2d Euclidean quantum gravity 
in the cases g) = (2,2m — 1), m = 2,3, . . ., can be described as double-scaling 
limits of one-matrix models with certain fine-tuned matrix potentials of order at 
least m + 1. Formally, the case m = 1, which corresponds to a somewhat degenerate 
(2, 1) conformal field theory with central charge c = —2, is then described by a 
special double-scaling limit of the purely Gaussian matrix model (see the review 
|34j). In this double-scaling limit one obtains for the so-called FZZT brane precisely 
the Airy function, see [351 EHl E] for recent discussions. 

Our model is not equivalent to 2d Euclidean topological quantum gravity, but 
is dual to it in a specific way. The arguments presented below suggest that our 
CDT string field theory can also be identified as a continuum theory associated 
with c = —2, but corresponding to the unconventional, "wrong" branch of the KPZ 
equation. - Recall that for a given conformal field theory coupled to Euclidean 2d 
quantum gravity, we have the KPZ formula 

for the susceptibility (for spherical topology). The parameter 7_ corresponds to 
the "right" choice of branch of the quadratic KPZ equation, i.e. the branch which 
leads to a weakly coupled Liouville theory as c ^ — oo. Choosing instead the other 
branch, one obtains 

-(l-c) + v/(l-c)(25-c) 7- .... 

^+ = 12 = -1^- ^''^ 

The interpretation of this 7_|. in terms of matrix models and geometry can be found 
in [38l |39| HOj , and for earlier related work see |^ W2\ . The simplest example is 
again given by c = —2: topological quantum gravity has 7_ = — 1 whose dual is 
the "wrong" 7+ = 1/2, which happens to be the value occurring generically in the 
theory of branched polymers (see, e.g., [lH for a discussion of why branched 
polymers and baby universes are generic and even dominant in many situations). 
While it is possible to describe 2d topological quantum gravity by a double-scaling 
limit of the Gaussian matrix model, the most natural geometric interpretation of 
the Gaussian matrix model is in terms of branched polymers, in the sense that the 
integral 



/ dM trM2" e" 



trM 



2 



/dM 



(37) 
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can be thought of as the gluing of a boundary of length n into a double-line branched 
polymer of length n. Since the branched polymers are also allowed to form closed 
loops, their partition function contains a sum over topologies 'en miniature', and one 
can indeed define a double-scaling limit of the model. When solving for the partition 
function in this limit, one obtains precisely our Z{g, A) of eq. f|T9|) n One should not 
jump to the conclusion that the two models are identical, since the CDT string field 
theory has a much richer structure of observables, with no obvious analogues in 
the branched-polymer set-up. Nevertheless, it is obvious that the partition function 
captures the essentials of the counting of branchings and joinings, which for the 
case of the CDT model is insensitive to the fact that the geometries are genuinely 
extended. (The latter is obvious from the nontrivial dependence on the boundary 
cosmo logical constants or, equivalently, the boundary lengths.) 

The fact that the CDT string field theory shares some properties of branched 
polymers is maybe less surprising in view of the fact that the original two-dimensional 
CDT model without branching can be mapped to a one-dimensional random-walk 
model [37j. The generalization implemented by the CDT string field theory cor- 
responds to adding "branches" to the random walks, resulting again in branched 
polymers. 

5 Discussion 

The central result of this paper is the derivation of the explicit, nonperturbative 
expressions ( |T9l) and ( l23l) for the partition function Z{g, A) and the Hartle-Hawking 
amplitude Wd{h, ...,/«) of the CDT string field theory, both incorporating the in- 
famous "sum over topologies" (2d spacetimes of all possible genera). It is rather 
remarkable that these sums can be performed and - with hindsight - in a man- 
ner which is technically not very involved. As was already the case for the two- 
dimensional CDT quantum gravity theory with fixed spacetime topology, the results 
are genuinely different from those of the corresponding purely Euclidean models. 
Nevertheless, as we have tried to argue above, from the point of view of conformal 
field theory, the CDT string field theory can probably be understood as the contin- 
uum theory "in the wrong branch" with conformal charge c = —2 and susceptibility 
7+ = 1/2, and thus as "dual" to topological quantum gravity in two dimensions. 

The nonperturbative aspects of our theory suffer from the same fundamental 
ambiguity as string theory, and which are rooted in the non-Borel summability of 
the perturbation series. Beyond the perturbative expansion, which unfortunately 
is only asymptotic, there is no real definition of the theory. We have managed to 
sum the asymptotic series and produce a closed formula for W{1), but like in non- 

^In the case of the branched polymer model, the parameters g and A appearing in the double- 
scaling limit are related to the "topology" of the branched polymer (i.e. to in the matrix 
model) and to the length of the polymer, respectively. We refer to [3S] for details on this work. 
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critical string theory with c > 0, i.e. unitary field theories coupled to 2d Euclidean 
quantum gravity, we lack a clear physical principle which would allow us to decide 
which nonperturbative completion of the perturbative expansion to choose. 

If one insists on a real partition function (in the Euclidean sector), it is natural 
to take F{z) = ttBiIz) in formulas f|T9|) and fl23|) . However, this choice is only 
unique within the matrix-model realization. While it is true that matrix models 
in non-critical string theory have been able to incorporate physics they were not 
originally designed to incorporate, like the physics of ZZ-hranes, we are not aware 
of any argument that would identify matrix models as the correct, nonperturbative 
definitions of continuum theories including a sum of genera, should they indeed 
exist. 

Our string-field theoretic model highlights in a particularly simple and trans- 
parent manner the limited amount of information contained in the perturbative 
expansion. When comparing explicitly the closed-form nonperturbative results (|T9l) 
and flM|) with their asymptotic expansions, one finds that the latter are only good 
approximations in a small range of their arguments t and /. This can simply be 
traced to the fact that the asymptotic expansion of for instance Bi(t~^/^), termi- 
nated after k terms, is only valid for t ^ A/k. In a similar vein, many aspects of 
the nonperturbative solution could not possibly have been guessed from the per- 
turbative series. Consider, for instance, the behaviour of the nonperturbative disk 
amplitude W{1) as a function of the boundary length / for fixed g and A. Each term 
in the perturbative expansion of W{1) falls off exponentially as e~^' for ^ 1 
and is positive. However, while the full nonperturbative function W{1) will ini- 
tially decrease with increasing /, as expected from each of the (positive) terms in its 
asymptotic expansion, it becomes oscillatory when / > l/{t\/\) = X/g. The same 
oscillatory behaviour occurs when t is increased while I and A are kept fixed, i.e. 
when the coupling g is increased. This oscillation is a genuinely nonperturbative 
effect, which is opposite to the behaviour of each term in the perturbative expansion 
of W{1) in powers of t and may be indicative of a phase transition to spacetimes 
completely dominated by topology changes. 
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